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Abstract 

The hybrid LU-SGS (Lower-Upper Symmetric 
Gauss-Seidel) algorithm has been added to a widely 
used series of 2D/3D Euler /Navier-Stokes solvers and 
has been demonstrated for a particular class of high- 
speed flows. A limited study has been conducted to 
compare the hybrid LU-SGS for approximate New- 
ton iteration and diagonalized Beam- Warming (DBW) 
schemes on a work and convergence history basis. This 
study showed that for the cases considered, the hybrid 
LU-SGS algorithm is more efficient and easier to imple- 
ment than the DBW scheme originally present in the 
code. The code has been validated for the hypersonic 
flow through two mutually perpendicular flat plates and 
then used to investigate the flow field in and around a 
simplified scramjet module gap seal configuration. Due 
to the similarities, the gap seal flow was compared to 
hypersonic corner flow at the same freestream condi- 
tions and Reynolds number. 

I. Introduction 

The hybrid LU-SGS (Lower-Upper Symmetric 
Gauss-Seidel) scheme has been added to the PARC 
2D/3D series of Euler/Navier-Stokes solvers. The 
PARC codes are based upon a central finite-difference 
scheme with fourth and second difference dissipative 
terms. The PARC codes originated as the AIR 2D/3D 
solvers by Pulliam (Ref. 1), which used the ADI al- 
gorithm presented by Beam and Warming (Ref. 2). 
To increase the computational efficiency and to reduce 
storage requirements, Pulliam (Ref. 3) diagonalized the 
block implicit operators of the original Beam-Warming 
ADI algorithm, resulting in the ARC 2D/3D codes. 
Cooper (Ref. 4) further modified the ARC codes to 
be used in a production environment for solving the 
propulsion oriented problems at the Arnold Engineer- 
ing Development Center, resulting in the PARC codes. 
The PARC codes are well validated and accepted and 
are in wide use by the aerodynamics community. 
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Through the use of a unique method of grid patching 
and the segmented application of boundary conditions, 
the PARC codes are very useful in solving wide varieties 
of flow problems. A full description of the PARC codes 
and their implementation of the diagonalized Beam- 
Warming algorithm can be found in Reference 4. 

In an attempt to further advance the practicality 
of these solvers, the present work investigates the ad- 
dition of a relatively new algorithm to the codes. An 
LU implicit scheme, proposed by Jameson and Turkel 
(Ref. 5) and improved upon and demonstrated for the 
Euler equations by Jameson and Yoon (Ref. 6) was 
used in conjunction with a multigrid strategy for the 
Euler equations by Jameson and Yoon (Ref. 7). From 
this scheme, the LU-SGS (also referred to as the LU- 
SSOR) scheme was developed and applied by Yoon and 
Jameson in References 8 and 9 for solutions to the Euler 
and Navier-Stokes equations. This scheme was shown 
to be very efficient and practical for the solution of 
the Navier-Stokes equations coupled in a fully implicit 
manner to the chemical transport equations by Shuen 
and Yoon (Ref. 10) in two dimensions and by Yu, Tsai 
and Shuen (Ref. 11) in three dimensions. Recently, 
this scheme was used by Park and Yoon (Ref. 12) to 
compute two dimensional flow in chemical and thermal 
non-equilibrium around a blunt body. This scheme has 
clear advantages when compared to others for implicitly 
solving fairly large systems of transport equations since 
it has the operational count and storage requirements 
of an explicit scheme, but still solves the equations im- 
plicitly and fully coupled. The following sections will 
show the derivation of the LU-SGS scheme, compare 
the LU-SGS and diagonalized Beam-Warming schemes 
in terms of convergence history on a work basis and use 
the LU-SGS scheme to analyze hypersonic corner and 
gap seal flowfields. 

EL Analysis 

The PARC2D/3D codes solve the Reynold’s av- 
eraged Navier-Stokes equations in a generalized curvi- 
linear coordinate system. The Navier-Stokes equations 
can be written in divergence form as 


the construction used in Reference 21, the approxi- 
mated jacobians are formed as 
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where the flux variables q and a generalized invis- 
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where 

ta,b,c = Pmax{\\ AiB ,c\) (8) 

where /3 > 1.0 and A A<B ,c are the eigenvalues of 
the flux jacobians. Therefore, using the conditioning in 
(7), the unfactored scheme may be written as 


where J is the jacobian of the coordinate transfor- 
mation and the contravariant velocities, u c , are defined 
in terms of the metric coefficients as 


u c = S x u + 6 y v + S z w (3) 

(The viscous fluxes are omitted here for brevity, 
and can be found in many References, such as Reference 
19.) Following the procedure outlined in Ref. 10 and 
augmented with results shown in References 20 and 21, 
the hybrid LU-SGS scheme can be derived from the 
following. A prototype implicit scheme for solving (1) 
in delta form may be written as 


[/ + At(D ( A + D V B + 0(0)^= -AtR (4) 

where D(, D v and D( are difference operators that 
approximate d{, d v and d( and A, B and C are the flux 
jacobians 


dE dF dG 

A = —,B = —,C = 

dq dq dq 


( 5 ) 


The PARC codes calculate the residual, R in (4), 
using central differences with added artificial dissipa- 
tion in the manner due to Jameson (Ref. 23). If the 
flux jacobians in (4) are approximated so that the re- 
sulting jacobians can be split into ”+” and matri- 
ces, where the eigenvalues of the ”+” matrices are non- 
negative and the eigenvalues of the matrices are 
non-positive, and then differenced according to their 
sign, (4) may be written as 


[I + At(V(A + -f A(A~ -f V v B + + A V B~ (6) 

+V ( C+ + A ( C-)] = -A tR 

The choice of conditioning used to form the ap- 
proximated jacobians is very important, and following 


[/(1-f At(r A + r B + r c ))-F At(A J+1 + 2? k+1 +C', +1 (9) 

-Af_i ~ B t-i ~ = -A tR 

where (j,k,l) refers to the grid point index. (When 
an index is omitted it is assumed to be at j,k or 1). The 
hybrid LU-SGS scheme is then found by factoring (9) 
into L,D and U operators so that 


LD~ 1 UAq— -AtR 

(10) 
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All of the calculations performed here are only con- 
cerned with obtaining the steady state flowfield, so an 
approximate Newton iteration formulation is formed by 
dividing (10) by A t and taking the limit as At ap- 
proaches infinity, resulting in 

LUAq — -(r A + r B + r c )R (14) 

L = [I(r A + r B + r c ) - A+_ t - - C+ ,] (15) 

U = [I{r A + r B + r c ) + AJ +1 + B^ +1 + C," J (16) 

The solution of (14) is fully vectorizable on planes 
where j+k+l= constant and is efficient since it requires 
only scalar diagonal inversions. The approach is to first 
solve for an intermediate vector using (15) by sweep- 
ing in the direction of planes of increasing j-f k+1 and 
then solve for A q using (16) by sweeping in the oppo- 
site direction. The efficiency in solving the chemical 
transport equations fully coupled to the Navier-Stokes 
equations has been demonstrated in References 10,11 
and 12 by using the LU-SGS scheme. The only modifi- 
cations to the implicit operator are to account for the 
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larger flux jacobians and the inclusion of the species 
production/destruction terms implicitly in the first op- 
erator, (15). Another advantage of this scheme is that 
when the approximate Newton iteration form is used, 
no choice of time step or Courant number is needed. 
This is very useful since for most problems it is not 
known a priori what time step or Courant number is 
needed to gain sufficient convergence. 

III. LU-SGS and DBW Comparison 

Since the present work has replaced the original 
diagonalized Beam- Warming scheme (hereafter referred 
to as DBW) with this algorithm it was deemed prudent 
to make comparisons regarding the relative efficiencies 
of the two schemes. A limited study was made for in- 
viscid 2D flows, and is summarized in the following sec- 
tions. In all the comparisons, only the implicit operator 
was changed; the residual calculations were performed 
in the same manner with the same coefficients of artifi- 
cial dissipation for both schemes. All calculations were 
performed on a CRAY-2 using the cft77 compiler. 

Ill.a: Moo =3 Blunt Body 

The inviscid flow over the windward portion of a 
right circular cylinder at a free stream Mach number of 
3 was calculated using the LU-SGS and DBW schemes. 
A 101 X 101 evenly spaced grid was used, resulting in 
approximately 62 percent of the grid points being in the 
subsonic portion of the flow. Figure 1 shows the com- 
puted flowfield Mach number contours which clearly in- 
dicates the curved bow shock and the expansion along 
the body surface from the stagnation streamline. For 
the DBW calculations a spatially varying time step was 
used by holding the Courant number (CN) constant at 
all the grid points while / 3 in (8) was linearly varied from 
it’s maximum, P max , to 1.0 over the first 500 iterations 
for the LU-SGS calculations. Figure 2 compares the 
convergence histories of the LU-SGS and DBW calcu- 
lations on a work basis, where one unit of work is equal 
to the cpu time for one iteration of the DBW scheme, 
for a range of Courant numbers and (3 max values. The 
Figure clearly shows the increasing rates of convergence 
with increasing CN for the DBW scheme and the rel- 
ative insensitivity to /3 max of the LU-SGS scheme. At 
CN’s higher than 2, the oscillations present in the be- 
ginning stages of the iteration persisted over a longer 
range, retarding the convergence to less than that for 
the CN=2 calculation. This Figure shows that the two 
schemes converged at roughly the same rates, but only 
when the DBW was run at a proper Courant number. 


For this case, the LU-SGS scheme was roughly 10 per- 
cent faster in time per iteration than the DBW scheme. 

IILb: Transonic Nozzle 

A simplified transonic nozzle flowfield was calcu- 
lated using the LU-SGS and DBW schemes. The ge- 
ometry was constructed by matching two cubic poly- 
nomials at the throat region with the inlet and exit 
heights determined from 1-D theory to provide an inlet 
Mach number of 0.2 and an exit Mach number of 3.0. 
The total temperature and pressure were held constant 
at the inflow boundary and the static pressure was held 
constant at the exit. A 101 X 51 evenly spaced grid was 
used for this calculation. Figure 3 shows the computed 
Mach number contours and geometry of the nozzle, 
showing that the flow accelerates smoothly through the 
throat and then rapidly becomes supersonic, where an 
excessive ammount of turning causes a weak shock near 
the upper nozzle lip. As in the previous case, spatially 
varying time stepping was used for the DBW scheme 
by holding the Courant number constant at all the grid 
points, while 0 was varied linearly from a maximum to 
1.0 over 500 iterations. Figure 4 compares the LU-SGS 
and DBW calculations for a range of Courant numbers 
(CN) and values of 0 max . As before, increasing the CN 
improved the convergence rate of the DBW scheme, but 
a limit was reached where beyond a CN=3.6, the DBW 
scheme diverged. Also, the LU-SGS was fairly insen- 
sitive to the 0 variation. For this case the LU-SGS 
was approximately 20 percent faster in time per itera- 
tion than the DBW scheme. The difference in efficiency 
when compared to the first case is attributable to the 
difference in the vector lengths (i.e.- grid sizes). Again, 
both schemes performed similarily, but only when the 
proper Courant number was used for the DBW scheme. 

IH.c: Moo =2.5 Supersonic Inlet 

The inviscid flow through a generic supersonic in- 
let at a freestream Mach number of 2.5 was calculated 
using the two schemes. The computational grid was 
uniformly spaced in both coordinate directions for this 
200 (axially) by 75 (normally) grid point calculation. 
Figure 5 shows the computed pressure contours in the 
inlet, indicating the passing of the ramp shock out of 
the domain and the cowl generated shock attenuation 
and propagation into the inlet, terminating in a large 
sonic region. The DBW calculation used a spatially 
varying time step and the LU-SGS calculation used a 
0 that varied from 1.5 to 1.0 over the first 100 itera- 
tions. The comparison of the DBW (with a CN=5.0) 
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and the LU-SGS schemes on a work basis is shown in 
Figure 6. For this case, the LU-SGS scheme performed 
approximately 40 percent faster per iteration than the 
DBW scheme. The Figure shows that the two schemes 
converged roughly the same at the early stages of the 
iteration, but the LU-SGS scheme converged faster in 
the later stages. 


Ill.d: Moo =3.5 Supersonic Inlet 

This case is identical to the previous case, except 
the freestream Mach number was increased to M=:3.5. 
The DBW scheme uses the same time stepping and 
Courant number as before, and the LU-SGS uses the 
same (3 variation. Figure 7 shows the computed pres- 
sure contours in the inlet, showing that the ramp shock 
does not get ingested and the cowl shock is strong 
enough to to travel throughout the inlet. Figure 8 com- 
pares the two schemes convergence behavior for this 
case on a work basis. The Figure shows, as in the pre- 
vious case, that the two schemes converge similarily 
in the early stages of the iterations while the LU-SGS 
scheme converges faster in the latter stages. As in the 
previous case, the LU-SGS scheme performed approx- 
imately 40 percent faster per iteration than the DBW 
scheme. 


Ill.e: Summary of Convergence Study 

The results of this study comparing the conver- 
gence rates and efficiencies of the LU-SGS and DBW 
schemes show that the LU-SGS scheme is slightly more 
efficient on a work basis than the DBW scheme. Since 
the LU-SGS vectorization is performed on planes in 3D 
(lines in 2D) where j+k+l=constant the efficiency of 
the scheme is dependent upon the grid dimension since 
the vector lengths will be short in the plane corners 
and longer in the center of the planes. Although the 
convergence rates were nearly the same on an iteration 
basis, the DBW scheme required more user effort to 
find the best time step or Courant number, while for 
the cases considered here, the LU-SGS was fairly insen- 
sitive to "tuning” of the implicit dissipation. Since one 
often does not know what the ’’best” time stepping (or 
Courant number) is needed, this study indicates that 
the LU-SGS for approximate Newton iteration scheme 
may be easier to use. 


IV: Corner and Gap Seal Flow Field Study 
Effort is currently underway to investigate the heat 


transfer and aerodynamic loads in and around a scram- 
jet module gap seal configuration. The seal system 
must prevent high temperature and pressure gases from 
leaking through the gaps between the articulating en- 
gine panel walls and stationary splitter panels (Ref. 
13). This system is crucial to the design of hypersonic 
airbreathing propulsion systems and presents a signifi- 
cant challenge to structural designers to build a viable 
seal system. To assist the design effort, the present re- 
search complements the structural analysis underway 
in an attempt to better understand the fluid dynamic 
and heat transfer phenomena in and around the gap 
seal system. 


IV.a: Hypersonic Corner Validation 


To obtain some level of confidence in the mod- 
elling of the gap seal flow it is necessary to validate the 
3D Navier-Stokes solver with experimental data taken 
from a class of flows similar to the gap seal flow. Fig- 
ure 9 shows the simplified geometry of a typical gap 
seal, showing that the geometry is similar to two mutu- 
ally perpendicular flat plates (corner) aligned with the 
freestream flow with a "groove” situated at the plate 
intersections. Due to this similiarity, the 3D code has 
been validated against the experimental data of hyper- 
sonic corner flow taken by Cresci (Ref. 14). These tests 
were conducted at a Mach number of 11.8 over a range 
of Reynolds numbers of 0.15 X 10® to 0.5 X 10® per 
foot. The data presented in Ref. 14 is axial pressures 
and heat transfer rates plotted against the hypersonic 
interaction parameter 


Rei 


(17) 


for different lateral locations measured from the 
corner. In addition, surface pressures, heat transfer 
rates and skin friction data are available (Ref. 15) as a 
function of lateral distance from the corner at various 
axial locations. Due to the simple geometry of this case, 
other researchers have successfuly used this data to val- 
idate space marching (Ref. 16) and Navier-Stokes (Ref. 
17) calculations. This data not only tests the predic- 
tion of the mutual interaction of the two plate boundary 
layers upon each other but also the self induced lead- 
ing edge pressures caused by the boundary layer growth 
close to the leading edges. To test the predictions of the 
2D self induced leading edge shocks and heat transfer 
rates, the 2D and 3D codes were run and compared with 
the experimental data of Ref. 14 (taken at a lateral 
distance sufficiently far from the corner to ensure two 
dimensionality) and with the theoretical pressures de- 
duced by Bertram (Ref. 18). The freestream conditions 
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were set at the lower Reynolds number and stagnation 
condtions from the experiment (Moo=ll. 8, Re=0.15 X 
10 6 /foot and To =1700 Rankine) while the wall tem- 
perature was held fixed at y“- =8.6543. Figures 10 and 
11 show the predicted and measured surface pressures 
and heat transfer rates versus x with the experimental 
and theoretical data (Ref. 18). A 36 X 41 grid was used 
with the grid clustered near the wall using a Robert’s 
transformation (Ref. 19) with 0 = 1.01 and axially near 
the leading edge with 0 = 1.1. Five grid points were ex- 
tended ahead of the leading edge to aid in correctly cap- 
turing the leading edge shock strength. The 3D code 
used a similar grid with 11 XY planes and symettric 
boundary condtions applied at the lateral boundaries. 
As seen in the Figure, the wall pressures were slightly 
over predicted while the heat transfer rates were pre- 
dicted very well. 

The results from the 2D study were then used to 
guide the 3D corner flow computations. Two grid sizes 
were used in the study (36 X 41 X 41 and 36 X 61 X 
61) for plate lenghts of 8 inches and plate widths of 2.5 
inches at a Reynolds number of 0.15 X 10 6 . The grids 
were clustered normally using the Robert’s transforma- 
tion in both planes with 0 =1.01 and axially near the 
leading edge using 0 = 1.1 with five points extended 
ahead of the leading edge. The convergence histories 
for the both grids using the LU-SGS scheme are shown 
along with a coarse grid calculation using the DBW 
scheme in Figure 12. The DBW calculation was made 
at a constant Courant number of 2.5, and contrary to 
the 2D inviscid results, the LU-SGS was only 3 percent 
faster in time per iteration than the DBW scheme. This 
Figure shows that for this case the LU-SGS computa- 
tions seemed to damp the initial transients more quickly 
than the DBW and that the convergence rate at the lat- 
ter stages of the iterations was faster for the LU-SGS 
scheme. At the inflow plane the flux variables were held 
constant at their freestream values, constant tempera- 
ture no-slip wall conditions were applied on the plates 
and grid extrapolation of the primitive variables was 
used at the outflow and lateral planes. Figures 13.a,b,c 
and d show the predicted and measured wall heat trans- 
fer rates measured laterally from the corner at differ- 
ent axial locations. Figures 14. a through e compare 
the predicted and measured wall pressures at different 
lateral locations plotted versus the hypersonic interac- 
tion parameter. Figure 15 shows spanwise surface pres- 
sure comparisons while Figure 16 shows the spanwise 
skin friction coefficients. Figures 13 a,b and d show re- 
markably good comparisons between the predicted and 
measured heat transfer rates, peak levels and peak lo- 
cations, while Figure 13. c shows the correctly predicted 
level but incorrectly predicted peak location. The pres- 
sure comparisons show good agreement near the corner 


(Figures 14. a, b and c) but get worse at farther distances 
from the corner (Figures 14.d and e). The experimen- 
tal data indicates that the plate pressure returns to the 
2D pressure at a distance of approximately 1.5 inches 
from the corner at \ = 5.0, but the present calculations 
only show a very gradual lateral reduction in pressure 
from the peak corner values. This behavior is similar to 
the computed results in Ref. 17 in that the near corner 
pressures are computed well, but the spanwise pressure 
relief is underpredicted. 

IV.b: Gap Seal and Corner Comparisons 

In an attempt to characterize the gap seal flow, 
comparisons were made between corner and gap seal 
flows at similar Reynolds and Mach numbers. Com- 
putations were made using the LU-SGS scheme for a 
Mach number of 10 at an altitude of 100,000 feet of 
the flow fields through a corner and a gap seal config- 
uration. At this altitude the freestream pressure and 
temperature are 23.085 psf and 418.79 Rankine yield- 
ing a Reynolds number per foot of 1.0 X 10 6 . The 
length of the plates and gap seal configurations were 
taken to be 5 inches with widths of 1 inch, resulting in 
a Reynolds number based on plate length of 420,000. 
The gap seal ”groove” (shown in Figure 9) was taken 
to be 1/8 by 1/8 inches and is considered to run along 
the entire length of the plates. The wall temperatures 
were held constant at 1200 Fahrenheit. The corner flow 
computation used a36X51X51 grid with the same 
clustering as the previous corner case with 5 points ex- 
tended ahead of the leading edge to help capture the 
leading edge shock strengths. The gap seal computa- 
tion used a rectangular grid with clustering centered 
along the gap seal faces and toward the plate walls as 
well as axial clustering near the leading edge with 5 
leading edge extension points. This computation used 
a base grid of 36 X 61 X 61 of which approximately 75 
percent of the grid points were active in the calculation. 
Figure 17 compares the convergence rates of the corner 
and gap seal computations using the LU-SGS scheme, 
showing that both cases converged monotonically and 
rapidly towards steady state. 

Figures 18.a,b and c compare the wall heat trans- 
fer coefficients along the gap seal floor and lateral wall 
with the corner only heat transfer rates for different 
axial locations. Examination of the Figures shows the 
peak heat transfer to the gap seal floor is reduced by 
approximately 4 times while the lateral wall heat trans- 
fer levels are approximately the same as the corner only 
levels. This redaction in heat transfer in the gap seal 
is due the the gap region filling up with low momen- 
tum and energy fluid caused by the viscous effects in 
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the gap. Figure 19 shows axial velocity contours near 
the end of the seal configuration in a plane perpendic- 
ular to the free stream velocity vector. This Figure 
shows the gap region filled with low velocity fluid and 
the existence of a vortical flow situated a short distance 
laterally from the gap seal. This vortical flow is caused 
by the rapid expelling of the flow from the gap seal 
to the outer lateral plate. Future computations will 
investigate whether this vortex will act as a barrier to 
keep the energetic outer flow from entering the gap seal 
region. 

V: Conclusions 

A hybrid LU-SGS scheme has been successfully im- 
plemented in the PARC2D/3D Navier-Stokes solvers 
and limited comparisons to the original algorithm in 
the code (diagonalized Beam- Warming ADI) have been 
made. The convergence comparisons indicated that the 
hybrid LU-SGS scheme for approximate Newton itera- 
tion is slightly more efficient on a work basis than the 
DBW scheme and requires less user effort to achieve 
an acceptable convergence rate. This is most likely due 
to a combination of the inherent implicit dissipation of 
the LU-SGS scheme and the approximate Newton iter- 
ation formulation, so that no time step specification is 
required. The scheme was demonstrated by validating 
the code against a =11.8 corner flow experiment 
and has been used to investigate the similarities be- 
tween gap seal and corner flow geometries. The valida- 
tion effort indicated that the code predicted the corner 
case heat transfer rates, peaks and peak locations very 
well while it predicted the wall pressures very well near 
the corner but poorly away from the corner. The gap 
seal/corner flow field comparisons indicated that the 
heat transfer rates to the floor of the gap seal were sig- 
nifigantly reduced when compared to the correspond- 
ing corner only flow and that the lateral heat trans- 
fer peaks were displaced laterally by approximately the 
gap width distance. Future work will investigate the 
gap seal/corner flows in the limit of large and small 
gap aspect ratios as well as the inclusion of inert gas 
’’blow by” cooling through the gap seal floor. 
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Figure 5 M x = 2.5 Supersonic Inlet: Pressnre Contours 


M = 2.5 Supersonic Inlet 

















Heat Transfer: x =7 ° Heat Transfer: *=4.6 




Figure 13 (b)Lateral Variation of Heat Transfer: x — 7-0 Figure 13 (d)Lateral Variation of Heat Transfer: x — 4.6 


Heat Transfer: ^=5.7 



Pressure: Z = 0.125 Inches 



Figure 13 (c)Lateral Variation of Heat Transfer: x — 5.7 Figure 14 (a)Aiial Variation of Wall Pressure: Z=0.125 

inches 
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Pressure: Z = 0.25 Inches 


Pressure: Z = 0.500 Inches 
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Figure 14 (d)Axial Variation of Wall Pressure: Z=0.500 

Figure 14 (b) Axial Variation of Wall Pressure: Z=0.250 inches 

inches 


Pressure: 7 = 0.375 Inches 


Pressure: Z= 1.000 Inches 
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Figure 14 (c)Ajrial Variation of Wall Pressure: Z=0.375 Figu rc 14 ( e ) Axial Variation of Wall Pressure: Z=1.0DO 
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Spanwise Pressure: x=5.0 Hypersonic Corner/Gap Seal 




Figure 15 Lateral Pressure Distribution at x = 5 0 Figure 17 Convergence History: Corner and Gap Seal Com- 

parison 


Skin Friction: x=5.0 



Figure 16 Lateral Skin Friction Distribution at x = 5.0 


Heat Transfer: Y = 0,X = 2 inches 



Figure 18 (a)Lateral Variation of Heat Transfer Rates, Cor- 
ner/Gap seal: X=2.0 inches 
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Figure 18 (b)Lateral Variation of Heat Transfer Rates, Cor- 
ner/Gap seal: X=4.0 inches 



Figure 18 (c)Lateral Variation of Heat Transfer Rates, Cor- 
ner/Gap seal: X=5.0 inches 
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